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1. Introduction 



L . 1.1. About these lecture notes. For centuries physics was a potent source pro- 

rK ' viding mathematics with interesting ideas and problems. In the last decades some- 

thing new started to happen: physicists started to provide mathematicians also 

i-Q ' with technical tools, methods, and solutions. This process seem to be especially 

strong in geometry and low-dimensional topology. It is enough to mention the 
mirror conjecture, Seiberg-Witten invariants, quantum knot invariants, etc. 

Mathematicians, however, en masse failed to learn modern physics. There seem 
to be two main obstructions. Firstly, there are few textbooks in modern physics 
written in terms accessible for mathematicians. Mathematicians and physicists 
speak two different languages, and a good "physical-mathematical dictionary" is 

l/"~j . missing 1 . Thus, to learn something from a physical textbook, a mathematician 

should start from a hard and time-consuming process of learning the physical jargon. 

^D , Secondly, mathematicians consider (and often rightly so) many physical methods 

and results to be non-rigorous and do not consider them seriously. In particular, 

f- ^ . path integrals still remain quite problematic from a mathematical point of view 

(due to some usually unclear measure aspects) , so mathematicians are reluctant to 
accept any results obtained by using path integrals. Yet, this technique may be put 
to good use, if at least as a tool to guess an answer to a mathematical problem. 
{""] In these notes I will focus on perturbative expansions of path integrals near a 

critical point of the action. This can be done by a standard physical technique of 
Feynman diagrams expansion, which is a useful book-keeping device for keeping 
track of all terms in such perturbative series. I will give a rigorous mathematical 
treatment of this technique in a finite dimensional case (when it actually belongs 
more to a course of multivariablc calculus than to physics), and then use a simple 
"dictionary" to translate these results to a general infinite dimensional case. 

As a result, we will obtain a recipe how to write Feynman diagram expansions 
for various physical theories. While in general an input of such a recipe includes 
path integrals, and thus is not well-defined mathematically, it may be used purely 
formally for producing Feynman diagram series with certain expected properties. 
A usual trick is then to "sweep under the carpet" all references to the underlying 
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2 MICHAEL POLYAK 

physical theory, keeping only the resulting series. Their expected properties often 
can be proved rigorously, directly from their definition. 

I will illustrate these ideas on the interesting example of the Chern-Simons the- 
ory, which leads to universal finite type invariants of knots and 3-manifolds. 

A word of caution: during the whole treatment I will brush aside all questions 
of measures, convergence, and such; see the discussion in Section 14.51 

1.2. Basics of classical and quantum field theories. The remaining part of 
this section is a brief sketch — on the physical level of rigor — of some basic notions 
and physical jargon used in the quantum field theory (QFT). Its purpose is to give 
a basic mathematical dictionary of QFT's and a motivation for our consideration 
of Gaussian-type integrals in this note. An impatient reader may skip it without 
much harm and pass directly to Section [21 Good introductions to field theories 
can be found e.g. in J2L EH : mathematical overview can be found in [5j; various 
topological aspects of QFT are well-presented in [23] ■ Very roughly, by a field theory 
one usually means the following. 

Given a space-time manifold X, one considers a space T of fields, which are 
functions of some kind on X (or, more generally, sections of bundles on X). A 
Lagrangian L : T — > R on T gives rise to the action functional S : T — > R defined 

by 



S(<j>) = / L(<j))dx. 
Jx 

In classical field theory one studies critical points of the action S ("classical 
trajectories of particles"). These fields can be found from the variation principle 
5S = 0, which is simply an infinite-dimensional version of a standard method for 
finding the critical points of a smooth function / : R — > R by solving M- = 0. 

In the quantum field theory one considers instead a partition function given by 
a path integral 

(1) Z = I e tksW V<t>. 

over the space of fields, for a constant k £ R and some formal measure T><f> on T . 
This is the point where mathematicians usually stop, since usually such measures 
are ill-defined. But let this not disturb us. 

In the quasi- classical limit k — > oo, the stationary phase method (see e.g. [S] 
and also Exercise 12. 5[) states that under some reasonable assumptions about the 
behavior of S this fast-oscillating integral localizes on the critical points of S, so 
one recovers the classical case. 

The expectation value (/} of an observable f : T — > R is 

For a collection f\ , . . . , f m of observables their correlation function is 

1 r n 

By solving a theory one usually means a calculation of these integrals or their 
asymptotics at k — > oo. 
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Increasingly often, due to a simpler behavior and better convergence properties, 
one considers instead the Euclidean partition function, equally well encoding phys- 
ical information (and related to P by a certain analytic continuation in the time 
domain, called Euclidean, or Wick, rotation): 

(2) Z= J e- kSW V(t). 

Since at present a general mathematical treatment of path integrals is lacking, 
we will first consider a finite dimensional case. 

1.3. Finite-dimensional version of QFT. Let us take T = M. d as the space of 
fields. An action S and observables fa are then just functions R d — ► E. For a 
constant k € R, consider the partition function Z = J Rd dxe~ ks ^ and the correla- 
tion functions (/i, . . . , f m ) = Z^ 1 j Rd dx e ~ kS ^ \\ i fi(x). We are interested in the 
behavior of Z and (f±, . . . , f m ) in the "quasi-classical limit" k — > oo. 

A well-known stationary phase method states that for large k the main contri- 
bution to Z and (f±, . . . , f m ) comes from some small neighborhoods of the points x 
where dS/dx — 0. Thus it suffices to study a behavior of Z and (/i, . . . , f m ) near 
such a point xq . Considering the Taylor expansions of S and /* in Xq (and noticing 
that the linear terms in the expansion of S vanish), after an appropriate changes 
of coordinates we arrive to the following problem: study integrals 

dx e-JM+WWpfj.) 

for some bilinear form A, higher order terms U(x), and monomials P{x) in the 
coordinates x l . 

Further in these notes we will calculate such integrals explicitly. To keep track 
of all terms appearing in these calculations, we will use Feynman diagrams as a 
simple book-keeping device. See the notes of Kazhdan in 9 for a more in-depth 
treatment. 

2. Finite-dimensional Feynman diagrams 

2.1. Gauss integrals. Recall a well-known formula for the Gauss integral (ob- 
tained by calculating the square of this integral in polar coordinates): 



Proposition 2.1. 



X dxe-^ = J*. 
-^ V a 



More generally, let A = (Ay) be a real d x d positive-definite matrix, x = 
(x 1 , . . . , x d ) the Euclidean coordinates in V = R d , and ( , ) : (R d )* xl^l the 
standard pairing (xi,x J ) — 8?. Then 

Proposition 2.2. 

(3) Z = f dx e~* <**•»> = fdet^- 

Indeed, by an orthogonal transformation (which does not change the integral) 
we can diagonalize A and apply the previous formula in each coordinate. 
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Remark 2.3. In a more formal setting, this may be considered as an equality for a 
positive-definite symmetric operator A : V — > V* from a d-dimensional vector space 
V to its dual (and ( , ) : V* x V -> E). Indeed, A induces det 4 : A d V -> A d V*, 
so that det A £ {k d V*)® 2 and (det A) 5 E \A d V\. Hence equality © with R d 
changed to V still makes sense if we consider both sides as elements of |A d y|. In a 
similar way, for C-valued symmetric operator A : V — > V* with a positive-definite 
ImA one has 

( dx e^ {Ax > x) = f det — 

where now both sides belong to |A d V|c- 

A more general form of equation J3J is obtained by adding a linear term — (6, a;) 
with b e (R d )* to the exponent: define Z\> by 



(4) Z b = dx g-i<Ar.») + <M). 

Then, by a change x — > .t — A -1 6 of coordinates, we obtain 
Proposition 2.4. 

(5) Z 6 = (det £\ ' ei^ A '^ = Z^ h ^ h \ 

Exercise 2.5. Verify the stationary phase method in the simplest case: use an 
appropriate change of coordinates to pass 

rP ax2 rP' ^ 

from / dx e k( ~^ +bx) to / dx e~^~ . 



What happens to a small e-neighborhood of the critical point xq = b/a under this 
change of coordinates? Conclude that in the limit k — > oo integration over a small 
neighborhood of Xo = b/a gives the same leading term in the expansion of this 
integral in powers of k, as integration over the whole of M. 

2.2. Correlation functions. The correlators (A, ■■•,/,«) of m functions /, : 
M. d — > R (also called m-point functions) are defined by plugging the product of 
these functions in the integrand and normalizing: 

(6) (fij2,...J m ) = ] r [dxe-^ A ^f l (x)... f m (x) . 

They may be computed using Z\>. Indeed, notice that 

^- f dx e -i<^>+<M> = J dx e -i<Ax,x)+(M)^ j 

hence for correlators of any (not necessary distinct) coordinate functions we have 
(I ) \x , . . . , x / — u n . . . u lm zji, i b — u n ... u lm e | b= „ 

where we denoted c\ = d/dbi. 

-Tr\mc! aro mi ran mr r rici l-l acci o n m o fviv . . . , . , . . . 



In particular, 2-point functions are given by the Hessian matrix M^{Z} ) JZq) 



with the matrix elements 

(8) {x\xi) = d i d j e^ A - Xh \_={A- 1 ) i 



1 6=0 ~ ^ ^' 
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Thus the bilinear pairing Sym (V*) — > R given by 2-point functions is just the 
pairing determined by A^ 1 . This explains the similarity of our notations for the 
2-point functions and ( , ) :V* xV ^R. 

For polynomials, or more generally, formal power series f±, . . . , f m in the coordi- 
nates we may apply (J7J) (with i\ = i% = • • • = i n = i for each monomial (x l ) n ) and 
then put the series back together, noting that each x % should be substituted by c^. 
This yields: 

Proposition 2.6. 

< 9 > <*■* «-a(s)-/-(b)'* <M "" , w 

2.3. Wick's theorem. Denote by A 1 ^ the matrix elements (A~ 1 )ij of A -1 . The 
key ingredient of the Feynman diagrams technique is Wick's theorem (see e.g. |22|1 
which we state in its simplest form: 



Theorem 2.7 (Wick). 



(10) d i± ...d im e~* 



0, m = 2n + 1 



where the sum is over all partitions (ji,J2),---, {jm—iijm) ^ n pairs of the set 
ii,i 2 ,- ■ ■ ,i m of indices. 

Proof. For each k, the expression di x . . .di k e^^ b ' A b ', considered as a function of 
b, is always of the form Pi 1 ...i k {b)e*( b ' A b ', where Pi 1 ...i h {b) is a polynomial. Each 
new derivative dj acts either on the polynomial part, or on the exponent, by the 
rule 

a,-(P(&)e4<*' A ~ l6 >) = dj(P(b))e^ b ' A ~ lb > +P(b)(J2A j %)e^ b < A ~ lb \ 

i 

so the polynomial part ^i...i m (fr) may be defined recursively by P${b) — 1 and 
(11) P n ... lm (b) = (d n +^A lll 6 l )P 2 ... 4m (6) = ••• 

i 

i i 

where 1 is the function identically equal to 1. We are interested in the constant 
term Pi 1 ...i m (0). Directly from (|ll|l we can make two observations. Firstly, if m 
is odd, P% 1 ...i m (b) contains only terms of odd degrees, in particular Pi 1 ..A m (0) = 0. 
Secondly, unless each derivative di k , k < m acts on the term ^2, i A lll bi in some 
Z-th, I > k, factor of the (|llf> . the evaluation at b = would give zero. Each such 
pair (ik,ii) contributes a factor of A lllk to the constant term of Pi 1 ...i m . These 
observations prove the theorem. □ 

It is convenient to extend (|10|) by linearity to arbitrary linear functions of the 
coordinates, note that in this case we may define the 2-point functions (/, g) by 
(f,A~ 1 g) in view of JSJ), and finally combine it with (JJJ into the following version 
of Wick's theorem: 
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Theorem 2.8 (Wick). Let fi(x), . . . , f m (x) be arbitrary linear functions of the 
coordinates Xi. Then all m-point functions vanish for odd m. For m — 2n one has 

(12) </l, ■ • -Jrn) = £</■„/«) • ■ ■ (/w-n/O, 

where the sum is over all pairings (i%, ii),. ■ ■ , (i m _i, i m ) ofl,...,m and the 2-point 
functions (fj,fk) are given by (fj,A~ l f k ). 

Remark 2.9. Another idea for a proof of Theorem 12.81 is the following. Note that 
both sides of Ijl2(l are symmetric functions of 1, . . . , m, so they may be considered 
as functions on 771-th symmetric power S m (V) of V = M. d . Thus it suffices to check 
(|12|l only for /i = ■ ■ • = f m = /; in this case it is obvious. 

Exercise 2.10. Check that the number of all pairings of l,...,2n is (2n)l/2 n nl. 
Calculate 

/>oo 

dx x m e x2/2 



using integration by parts and Proposition 12. II Calculate 

d" 



dx' 



. e *7» 



substituting x 2 /2 instead of x in the Taylor series expansion of e x . Compare these 
expressions and explain how are they related to the above number of pairings. 

Exercise 2.11. Find formulas for the 4-point functions (xi, x\, x 2 , X3) and (x\, X\, x\ 1 xz)- 

2.4. First Feynman graphs. It is convenient to represent each term 

\Jil : Hi) ■ ■ ■ \Ji,n-i 5 Jim) 

in Wick's formula IJ12JI by a simple graph. Indeed, consider m points, with the fc-th 
point representing ft- . A pairing of 1, . . . , 2n gives a natural way to connect these 
points by n edges, with an edge (a propagator in the physical jargon) e = (j, k) 
representing A^ 1 = (fj,A~ 1 f k ). Equation l(T2"|l becomes then 

(13) (/i,...,/ m >=£ n A e^ 

r eecdgos(r) 
where the sum is over all univalent graphs as above. 

Example 2.12. An application of equation l|13fl for n = 2 (see Figure^,) gives the 
following: 

{xx, x 2 , x 3 , Xi ) = A 12 A 34 + A 13 A 24 + A 14 A 23 , 

(x u x ll x 2 ,x 2 )=A 11 A 22 + 2A 12 A 12 1 

(a;i,xi, £1,2:1} = 3A n A n . 

2.5. Adding a potential. The above computations may be further generalized by 
adding a potential function U{x) (with some small parameter h = k^ 1 ) to (Ax,x) 
in the definition of Zq. Namely, define 2 Zjj by 

(14) Z v = f dxe-^ Ax ' x)+hu{x) . 

Applying © for / = e hu ^ we get: 

Again, let me remind that we ignore problems of convergence: for most U(x) this integral will 
be divergent! 
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1 3 1 3 1 ,3 



4 2 4 2 

a b 

Figure 1. Terms of {x\, X2,x^,X4) and graphs of degree two 

Proposition 2.13. 

(15) Z u = Z Q e hU ^e^ A - X \ =Q . 
Correlation functions (/i, . . . , f m )u are defined similarly to JBJ: 

(16) (A, / a , . . . J m ) v = ±- J dx e -i(A^)+m( x)fi{x) fk{x y 

Using © once again, we get 
Proposition 2.14. 

(IT) (/„*,...,/„)„ = fe-W A (»).../. (|) .»**-•• 

2.6. A cubic potential. Consider the important example of a cubic potential 
function U(x) = '^2,UijkX t x 3 x k . Let us compute the expansion of the partition 
function 114|) in power series in h. The coefficient of h n in the expansion of (|15|l is 



6=0 






6=0 



Let us start with the lowest degrees. By Wick's theorem, the coefficient of h vanishes 
and the coefficient of h 2 is given by 

( 18 ) frE E u ijk u i , j , v a i d J d k a i ,d j ,d k ,eW> A - lb '>\ b=0 = 



§E E u l3k ue rk >Y. A%ll2Al3liA ^^ 

i,j,k i' ,j' ,k' 

where the last sum is over all pairings (ii,i2),- ■ ■ ,(i5,ie) of i,j,k,i',j',k'. We 
may again encode these pairings by labelled graphs, connecting 6 vertices labelled 
by i,j,k,i',j',k' by three edges («i,i2),(*3)U),(«5,i6) representing A ni2 A 1314 - A 1516 . 
This time, however, we have an additional factor Vi^kUi'j'W- To represent Uijk 
graphically, let us glue the triple (i,j, k) of univalent vertices in a trivalent vertex; 
to preserve the labels, we can write them on the ends of the edges meeting in this 
new vertex (i.e., on the star of the vertex). Similarly, we represent Ui>j'k> by gluing 
the remaining triple (i',j',k') of univalent vertices into a second trivalent vertex. 
Thus for each of the 6!/(2 3 3!) = 15 pairings of i,j,k,i',j',k' we end up with a 
graph with two trivalent vertices; we get 6 copies of the 0-graph and 9 copies of 
the dumbbell graph shown in Figure \T]p. 

Note, however, that each of these labelled graphs is considered up to its auto- 
morphisms, i.e. maps of a graph onto itself, mapping edges to edges and vertices 
to vertices and preserving the incidence relation. Indeed, while the application of 
an automorphism changes the labels, it preserves their pairing (edges) and the way 
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Figure 2. Degree two graphs with two legs 

they are united in triples (vertices), thus corresponds to the same term in the right 
hand side of lfT%|) . Instead of summing over the automorphism classes of graphs, 
we may sum over all labelled graphs, but divide the term corresponding to a graph 
r by the number | AutT| of its automorphisms. E.g., for the O-graph of Figure^) 
| AutT| = 12, and twelve copies of this graph (which differ only by transpositions 
of the labels) all give the same terms U ijkU i> j' w A %% A" A kk . 

Also, when summing the resulting expressions over all indices, note that the 
terms corresponding to i, j, k, i' ', j' , kl and to i', j', k',i,j,k are the same (which will 
cancel out with 1/2! in front of the sum). Hence, we may write the coefficient of 
h 2 in the following form: 



*£n^rr e II MI* 1 . 



lAutri 

r labels v e 

Here the sum is over all trivalent graphs with two vertices and labellings of their 
edges, U v — Vijk for a vertex v with the labels i,j,k of the adjacent edges, and 
A^ 1 — A^ for an edge e with labels i,j. 

Exercise 2.15. Calculate the number of automorphisms of the dumbbell graph of 
Figure HJ). 

In general, for the coefficient of h n we get the same formula, but with the sum- 
mation being over all labelled trivalent graphs with n vertices. 

2.7. Correlators for a cubic potential. We may treat m-point functions in a 
similar way. Let us first consider the power series expansion in h of Zjj (x 11 , . . . , x lm ) u . 
The coefficient of h n is 



^ EMM] d h ...d im e^M-6) 

\i,j,k 



6=0 



Thus it may again be presented by a sum over labelled graphs, with the only 
difference being that now in addition to n trivalent vertices these graphs also have 
m ordered legs (i.e. univalent vertices) labelled by ix,...,i m . See Figure for 
graphs representing the coefficient of H 2 in Zu(x 1 ,x 2 )u. 

However, not all of these graphs will enter in the expression for (a;* 1 , . . . , x lm )u, 
since we should now divide this sum over graphs by Z\j (represented by a similar 
sum, but over graphs with no legs). This will remove all vacuum diagrams, i.e. all 
graphs which contain some component with no legs. Indeed, the term corresponding 
to a non-connected graph is a product of terms corresponding to each connected 
component. Each component with no legs appears also in the expansion of Z\j and 
thus will cancel out after we divide by Zjj . For example, the first graph of Figure 
12 contains a vacuum 0-graph component. But it also appears in the expansion of 
Z\j (see Figure ^p) . Thus the corresponding factor cancels out after division by 
Zu- The same happens with the second graph of Figure As a result, only the 
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last two graphs of Figure El will contribute to the coefficient of ti 2 in the expansion 
of {x l ,x 2 ) v . 

Example 2.16 ("A finite dimensional 3 -theory"). Take E/y& = SijSjk, i.e. U — 
'}2, i {x l Y . Note that since all ends of edges meeting in a vertex are labelled by the 
same index, we may instead label the vertices. Thus the calculation rules are quite 
simple: we count uni-trivalent graphs; a vertex represents a sum Ei over its labels; 
an edge with the ends labelled by i,j represents A % K The coefficients of h 2 in Zjj 
and in (x 1 1 x 2 )jj are given by 



Z £ Q(A lJ ) 3 + 9A ij A u A jj , 

' 9A u A 2j A u A» + 6A u A^(A^f 



i,j 



i,0 



respectively. We can identify these terms with two graphs of Figure ^> and two 
last graphs of Figure |3 respectively. The first two graphs of Figure |21 represent 



i,0 



6A 12 (A ij ) 3 + §A 12 A l3 A il A" 



These terms do appear in the h 2 coefficient of Zjj(x 1 1 x 2 ), but cancel out after we 
divide it by Z v =Z (l + h 2 Ey ( 6 (^ y ) 3 + 9A^A U A») + ...). 

2.8. General Feynman graphs. It is now clear how to generalize the above re- 
sults to the case of a general potential U(x): the k-th degree term Ui 1 ...i k (c n . . . x lk 
of U will lead to an appearance of fc-valent vertices representing factors Ui x ...i k ■ We 
will call such a vertex an internal vertex. We assume that there are no linear and 
quadratic terms in the potential, so further we will always assume that all internal 
vertices of any Feynman graph T are of valence > 3; denote their number by |T|. 
Denote by T° the set of all graphs with no legs. Also, for m > 1, denote by T m the 
set of all non-vacuum (i.e. such that each connected component has at least one 
leg) graphs with m ordered legs. 

Denoting U v = Ui 1 ..,i k for an internal vertex v with the labels i%, . . . ,ik of adja- 
cent edges, and A~ x — A %: > for an edge e with its ends labelled by i,j, we get 

Proposition 2.17. 

( 19 ) ** = *>£ j^EII^IK 1 - 

rer° ' ' labels v e 

Note that instead of performing the internal summation over all labellings, one 
may include the summation over labels of the star of a vertex into the weight of 
this vertex. 

In a similar way, for m-point functions we get 

Proposition 2.18. For even m, 

( 2 °) (^,...^)u= £ ^- £ II 17 ' IK 1 - 

Ter m ' ' labels v e 

where the sum is over all labelled graphs T with m legs labelled by i\, . . . , i m . 

Again, we may include the summation over the labels of the star of an internal 
vertex into the weight of this vertex. 
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2.9. Weights of graphs. Let us reformulate the above results using a general 
notion of weights of graphs. 

Let V be a vector space. A weight system is a collection (a, {uk}kL s ) of a G 
Sym 2 (V0 and u k G Sym fe (F*). A weight system W defines a weight W T : (V*)® m -» 
R of a graph T G T™ 1 in the following way. Assign u k G Sym (V*) to each internal 
vertex v of valence k, associating each copy of V* with (an end of) an edge. Also, 
to the i-th leg of T, i = 1, . . . , m assign some fi G V*. Now, for each edge contract 
two copies of V* associated to its ends using a G Sym 2 (V). After all copies of V* 
get contracted, we obtain a number Wr(/i, ■ ■ ■ , f m ) G K. 

In our case, a bilinear form A^ 1 and a potential KU{x) determine a weight 
system in an obvious way: set a = A^ 1 and let u v to be the degree k part of hU(x). 
These rules of computing the weights corresponding to a physical theory are called 
Feynman rules. 

Formulas 1)19(1 and l|2()[) above can be reformulated in these terms as 



%> = Z »E TT^rl^ 



rer° 

lAutr 



(21) ^ , 

(/l, ■ • ■ ,fm)u = 2^ I AntFl ^" 1 ^ 1 ' " " •'■^ T ")- 

rgr" 

Exercise 2.19 (Finite dimensional (/> 4 -theory) . Consider a potential C/ = X^f^ 1 ) 4 - 
Formulate the Feynman rules. Find the graphs which contribute to the coefficient 
of h 2 of Zjj and compute their coefficients. Do the same for (x 1 ,x 2 )u- Draw the 
graph representing ^ A 12 (A 11 ) 2 ; does it appear in the expansion of (x 1 ,x 2 )u and 
why? 

2.10. Free energy: taking the logarithm. The summation in equation l|21|l is 
over all graphs in F°, which are plenty. Denote by r" on „ the subset of all connected 
graphs in T°. There is a simple way to leave only a sum over graphs in r° onrl , 
namely to take the logarithm of the partition function (called the free energy in the 
physical literature): 

Proposition 2.20. Let W be a weight system. Then 




Autri / v lAutr 

Proof. Let us compare the terms of the power series expansion for the right hand 
side with the terms in the left hand side: 

where the sum is over all fc, n i: and distinct Ti G r° OTm , i = l,...,k. Consider 
T = (Ti) ni . . . (T k ) nk G r°. Since in addition to automorphisms of each r, there 
are also automorphisms of T interchanging the n, copies of Tj, we have | Aut T| = 
m ! ... rife! | Aut Ti | . . . | Aut T k \ ■ Also, any weight system satisfies Wtt" = Wr> Wr» , 
hence W T = W^ . . . W^ . The proposition follows. □ 

Exercise 2.21. Formulate and prove a similar statement for graphs with legs. 
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Remark 2.22. It is possible to restrict the class of graphs to 1-connected (in the 
physical literature usually called 1-point irreducible, or 1PI for short) graphs. A 
graph is 1-connected, if it remains connected after a removal of any one of its edges. 
This involves a passage to a so-called effective action, which I will not discuss here in 
details. Mathematically, it simply means an application of a Legendrian transform 
(a discrete version of a Fourier transform): if z(b) = log(Z(,) is given by the sum 
over all connected graphs as in Proposition 12 .201 then z(x) = (b,x) — zip) is given 
by a similar sum over all 1PI graphs (and b(x) may be recovered as dz/dx). 

3. Gauge theories and gauge fixing 

3.1. Gauge fixing. All calculations of the previous section dealt only with the case 
of a non-degenerate bilinear form A; in particular, the critical points of the action 
S(x) had to be isolated (see Section IT"3"|) . However, gauge theories present a large 
class of examples when it is not so. Suppose that we have an /-dimensional group 
of symmetries, i.e. the Lagrangian is invariant under a (free, proper, isometric) 
action of an /-dimensional Lie group G. Then instead of isolated critical points we 
have critical orbits, so A has / degenerate directions and the technique of Gauss 
integration can not be applied. 

Let us try to calculate the partition and correlation functions without a superflu- 
ous integration over the orbits of G. In other words, we wish to reduce integrals of 
G-invariant functions on X to integrals on the quotient space X — X/G of G-orbits. 
For this purpose, starting from a G-invariant measure on X we should desintegrate 
it as the Haar measure on the orbits over some "quotient measure" jj, on the base 
X. 

If G is compact then jj, is the standard push-forward of \i. For example, if / is a 
rotationally invariant function on K 2 , we can take the pair of polar coordinates (r, (j>) 
as coordinates in the quotient space X and the orbit, respectively. The measure on 
X in this case is 2irr dr and we get the following elementary formula: 

/>OC 

/(|x|) d 2 x = 2ir f(r)r dr. 

Jo 

If G is a locally compact group acting properly on X then /} can be defined by 
the property 



fx(Y) = / \YnG(x)\dji(x), 
J x 

where G(x) is the fiber over x G X and | • | is the Haar measure on it. In this case the 
integral J x f dx in question is infinite, but it can be formally defined ("regularized") 
as Jy fdx. 

A standard physical procedure for the desintegration that can be applied also 
to a non-locally compact gauge group is called a gauge fixing (see e.g., J22]); it 
goes as follows. Suppose that / : X — > R is G-invariant, i.e. f(gx) = f(x) for 
all x G X, g G G. Choose a (local) section s : X — > X which intersects each 
orbit of G exactly once. Suppose that it is defined by I independent equations 
F x (x) = • •• = F l (x) = for some F : X — > M. . Firstly, we want to count each G- 
orbit only once. This is simple to arrange by inserting an /-dimensional (5-function 
S l (F(x)) in the integrand. Secondly, we want to take into account the volume of a 
G-orbit passing through x, so we should count each orbit with a certain Jacobian 
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factor J(x) (called the Faddeev- Popov determinant) . How should one define such a 
factor? We wish to have 

f(x) dx= I f(x)J(x)6 l (F(x)) dx. 
x Jx 

Rewriting the right hand side to include an additional integration over G and 
noticing that both f(x) and J{x) are G-invariant, we get 

f{x) dx= I f(x)J{x)S l {F(x) dx = 

i A JX 

dx f dg f(x)J(x)S l (F(gx)) = f dx f(x)J(x) f dg 5 l {F{gx)). 

IX JG JX JG 

Thus we see that we should define J(x) by 



J(x) / dgS l (F(gx)) = l, 

JG 

where dg is the left G-invariant measure on G. Thus, the Faddeev-Popov deter- 
minant plays the role of Jacobian for a change of coordinates from x to (s(x),g). 
Example in H3. 31 below provides a good illustration. 

Remark 3.1. A formal coordinate- free way to define J(x) is as follows. The section 
s : X — > X determines a push-forward s* : T X X — > T X X of the tangent spaces. The 
tangent space T X X thus decomposes as s* (Bi ■ T x X(Bg — * T X X, where i : g — > T X X 
is the tangent space to the orbit, generated by the Lie algebra g of G. The Jacobian 
J(x) may be then defined as J{x) = det(s* © i). 

Remark 3.2. Equivalently, one may note that the tangent space to the fiber at 
x 6 s may be identified with g, to directly set J(x) — det A, where A = (§&-) and 
{fl J }i = i is a set of generators of the Lie algebra g of G, see e.g. |3]. I.e., J(x) is the 
inverse ratio of the volume element of g and its image in R' under the action of G 
composed with F . 

Indeed, since F has a unique zero on each orbit and since (due to the presence 
of the delta- function) we integrate only near the section s, we can use F as a local 
coordinate in the fiber over x. Making a formal change of variables from g to F we 
get 

J(x)- 1 = J dg S\F{gx)) = J dF S l (F(gx)) det (^) = det ^ 

Calculating ( — q' ) __ at a point x € s and identifying the tangent space to the 
fiber with g, we obtain (ff )| F=0 = (ffr)- 

Exercise 3.3. Let us return to the simple example of a rotationally invariant function 
f{x\,X2) — /(M) on R 2 , using this time the gauge-fixing procedure. The group 
G = S 1 acts by rotations: <px = e^x and the (normalized) measure on G is ^-dtfi. 
We should use the positive xi-axis for a section s, so we may take e.g. F = X2- A 
slight complication is that the equation X2 = defines the whole a;i-axis and not 
only its positive half, so each fiber of G intersects it twice and not once. This can 
be taken care of, either by dividing the resulting gauge-fixed integral by two, or 
by restricting its domain of integration to the right half-plane R^ in R 2 . In any 



F=0 
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case, using x 2 instead of <f> as a local coordinate in the fiber Gx near i£swe get 
d<fi = d[ arctan(a;2/a;i)) = xi\x\~ 2 dx2 Thus for x £ s we have 

J(x)- 1 = I ' 5{F{4>x)) ^# = ±- [ 6(x 2 )xi\x\- 2 dx2 = J-, 

so J(|x|) = 27r|x| as expected and 

f(\x\)d 2 x= ( f{\x\)2ir\x\5{x 2 )dx 1 dx 2 = 2Tr [ f(r)rdr. 
m\ Jo 

3.2. Faddeev-Popov ghosts. After performing the gauge-fixing, we are left with 
the gauged-fixed partition function 

Z GF = [ dxe-^ Ax ' x) 6 l (F(x))detA. 



We would like to make it into an integral of the type we have been studying be- 
fore. We have two problems: to include 8(F(x)) det A in the exponent (i.e., in the 
Lagrangian) and — more importantly — to make A into a non-degenerate bilinear 
form. 

The (5-function is easy to write as an exponent using the Fourier transform: 

6 l (F{x)) = {2ir)- 1 f d^e l ^ F(x)) . 



The gauge variables £ (called Lagrange multipliers) supplement the variables x, and 
the quadratic part of (£, F(x)) supplements {Ax,x) so that the quadratic part Ap 
of the gauge-fixed Lagrangian is non-degenerate. 

The det A term is somewhat more complicated; it can be also represented as 
a Gaussian integral, but over anti- commuting variables c = (c 1 , . . . ,c l ) and c = 
(c 1 , . . . , c'), called Faddeev-Popov ghosts. Thus 

(22) cV + cV = cV + e'e 1 = c¥ + c 7 '? = 

There are standard rules of integration over anti-commuting variables (known to 
mathematicians as the Berezin integral, see e.g. J23I Chapter 33] and J13II16J ). The 
ones relevant for us are 

/ c l dc j = / ? dc j = 8^ and / 1 <k? = f 1 dc) = 0. 

The multiple integration (over e.g., dc = dc l . . .dc 1 ) is defined by iteration. One 
may show that this implies (see the Exercise below) that for any matrix A 

e <E < Ac> dcdc = det A. 

Exercise 3.4. Let I — 1 and define the exponent e Acc by the corresponding power 
series. Use the commutation relations (1221) to verify that only the two first terms 
of this expansion do not vanish. Now, use the integration rules to deduce that 
Je XEc dcdc = /(1 + Acc) dcdc = A. 

Thus we may rewrite Zqf by adding to the Lagrangian the gauge-fixing term 
and the ghost term: 



gf = I dxd^dede e 2 



±(Ax,x) + {c,Ac)+i(£,F(x)) 
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At this stage we may again apply the Feynman diagram expansion to the gauge- 
fixed Lagrangian. The Feynman rules change in an obvious fashion. The quadratic 
form now consists of two parts: Ap and A, so there are two types of edges. The first 
type presents Ap, with the labels x % and £ l at the ends. The second type presents 
A, with the labels c 1 and <? at the ends. Note that since A is not symmetric, these 
edges are directed. Also, there are new vertices, presenting all higher degree terms 
of the Lagrangian (in particular some where edges of both types meet). An example 
of the Chern-Simons theory will be provided in Section |SJ 



2 . In the coordinates (x±, X\, £2,22) on 
— > Xxi. Let us take A — 2121 as an 

' X-2 X-2 



3.3. An example of gauge-fixing. Let us illustrate the idea of gauge-fixing on 
an example of the standard C*-action on 1 
C 2 the gauge group acts by Xi — > Xxi, x 
invariant function. 

Of course, the orbit space CP 1 is quite simple and an appropriate measure 
on CP 1 is well known; in the coordinates z — Xi/x%, z = X\/x2 it is given by 
dzdzjil + zz) 2 We are thus interested in 

f dzdz 1 „= 

(23) Za ' = JoTW e • 

Let us pretend, however, that we do not know this and proceed with the gauge- fixing 
method instead. 

The invariant measure on C 2 is dx\dx\dx 2 dx2 / (x\X\ +X2X2) 2 '■ In a gauge F = 
we have 

dxidx2dxidx2 



Zgi 



(X1X1 +X2X2) 2 



^S 2 (F(x,x))detA, 



where A 



x\F Xl +x 2 Fx 2 

x\F Xl + x 2 F X2 ^n, 

E.g., for F = X2 — 1 we get 8 2 {\x2 — 1|) and det A = X2S2 



XlFx! 

x\F x 



X2Fjc 2 
X2F X2 



Exercise 3.5 (Different gauges give the same result). Consider F = x% — 1. Show 



that S 2 ( 



1 = 



2 S 2 (\x 2 — 1|) and det A = aa(x 2 X2) a - Check that the 



dependence on a in Zqp cancels out, thus gives the same result as F = x% — 1. 
Show that it coincides with formula (|23() . 

Finally, let us check that while the initial quadratic form A is degenerate, the 
supplemented quadratic form Ap is indeed non-degenerate. It is convenient to make 
a coordinate change x[ = xi, x' 2 = £2 — 1. Using a Fourier transform we get 



5 2 (x 2 - 1) = (2tt)- 



d£d£e 



i(£x 2 -£x' 2 ) 



Also, we have — 

X2 



c i 12^Li(~^) n x 2- We can now compute A and Ap; in the 



coordinates (x^ , x' x , x' 2 , x' 2 ) and (xi, Xi, x' 2 , x' 2 , £,,£,) , respectively, we have: 



A 



10 

10 







-4, 



1 

1 


i 
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4. Infinite dimensional case 

4.1. The dictionary. Path integrals are generally badly defined, so instead of try- 
ing to deduce the relevant results rigorously, we will just provide a basic dictionary 
to translate the finite dimensional results to the infinite dimensional case. 

The main change is that instead of the discrete set i G {1, . . . , d} of indices we 
now have a continuous variable x € M n (say, in K"), so we have to change all 
related notions accordingly. The sum over i becomes an integral over x. Vectors 
x = x(i) = (x 1 , . . . , x d ) and b — b(i) become fields cf> = <fi(x) and J(x). A quadratic 
form A = A(i,j) becomes an integral kernel K = K(x,y). Pairings {Ax,x) = 
Y^i j x l Aij xi and (b,x) = Yi^ %x% become (K(f),<f)) = J dxdy 4>{x)K{x,y)4>(y) and 
(J, <j)) — J dx J{x)4>(x) respectively. The partition function Z\, defined by (@J 
becomes a path integral Zj over the space T of fields 



Zj = / V<\> e 

The inverse A -1 of A defined by ^ fc AikA k; > = Sj corresponds now to the inverse 
G = K- 1 of K defined by 



dz K{x 1 z)G(z, y) = 6(x — y). 

Formula © for Z\> then translates into 

Zj = Z e*< J ' GJ >. 
Correlators (x H , . . . , x lm } defined by 10 become now m-point functions 

(<l>(xi), ..., <jy{x m )) = -i- f-Dtf, e-W+'+Uixi) ■ . ■ <p(x m ). 
^o J 

4.2. Functional derivation. A counterpart of the derivatives d/dx l is given by 
the functional derivatives 5/5<fi(x). The theory of functional derivation is well- 
presented in many places (see e.g. |l(Jp. so I will just briefly recall the main notions. 
Let F(<j>) be a functional. If the differential 

DFMW-UmW + rt-™ 

SF _ 



can be represented as J p(x)h(x)dx for some function h(x), then we define 
h(x). In general, the functional derivative 5F/S<p(x) is the distribution representing 
the differential of F at <\>. The reader can entertain himself by making sense of the 
following formulas, which show that its properties are similar to usual derivatives: 

-4>(y) = S(x-y), 



5<f>( 



x 



' ] (F(cb)H(4>)) = t^(F(4>)) ■ H{4>) + F{4>) ■ *(H{<I>)). 



S(f>(x) §4>(x) 5<fi{x) 

Example 4.1. 

5 e {M) - 5 e J dxJ(x)<j>(x) = My} e J dxJ(x)4>(x) _ 0( y ) e <J,0>. 



5J(y) 6J(y) 
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Exercise 4.2. Consider a (symmetric) potential function 



Y^ 



1 



(24) U{4>) = > — - I dxi . ..dx n U n (xi,. .. ,X n )cf>(xi) . . .</)(x n ) 



Prove that 
5 



■U((f>) = ^2~ j dxi.. .dx n U n+1 (y,x 1} .. .,x n )(j>(xi) ■ ..4>{x n ). 



H(y) 

The inverse G(x, y) can be written as a Hessian, similarly to equation (jSJ) for 

-1 c c 

Z dJ(a;) d J(y) J=0 
More generally, for m-point functions we have, similarly to @, 

1 5 5 



(<f>(x !),..., </){x m )) = 



.7=0 



Z <5J(xi) ' " 5J(x m ) 

4.3. Wick's theorem and Feynman graphs. Wick's theorem now states that, 
similarly to I jlOjl , 

T77 r • ■ • TT7 \ e2 = / _, ^{Xin x i2) • ■ ■ ^ Jr \Xi m _ 1 , Xi m J, 

dJ(xi) oJ(x m ) J=0 ^^ 

where the sum is over all pairings {i\, i%) . . . (i m -i, im) of 1, ... , m. Just as in the 
finite dimensional case, we may encode each pairing by a graph with m univalent 
vertices labelled by 1, . . . , m, and edges connecting vertices i\ with i%, . . . , and 
i-m-i with i m presenting the factors of G. 

Let us add a potential (|24|l to the action and define 

Z u= f V(j)e -m(K^4>)+hU(4>)). 

Then, similarly to l|15|). we have 

Z u = Z e hu ^ei< J ' GJ >\ J=0 . 

Using again the Wick's theorem, we can rewrite the latter expression in terms of 
Feynman graphs to get 

,25) *-?E*jLn".n«« 

where the integral is over all labellings of the ends of edges, U v = U(xi, . . . , Xk) for a 
fc-valent vertex with the labels X\, . . . , Xf. of the adjacent edges, and G e — G(xi,Xj) 
for an edge with labels Xi, Xj. Sometimes it is convenient to include the integration 
over the labels of the star of a vertex into the weight of this vertex. 

4.4. An example: 4 -theory. Let us write down the Feynman rules for a poten- 
tial U{4>) = J dx </> 4 (x). Firstly, the relevant graphs have vertices of valence one 
or four. Secondly, all edges adjacent to a vertex should be labelled by the same x, 
so we may instead label the vertices. An edge with labels x, y represents G{x, y) 
and (including the integration over the vertex labels into the weights of vertices) 
an relabelled vertex represents J dx. The linear term in the power series expansion 
of (#i, X2)u should correspond to non-vacuum graphs with two legs, labelled by X\ 
and X2, and one 4-valent vertex. There is only one such graph, see Figure El It 
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Figure 3. Graphs appearing in the </> 4 -theory 

represents J dx G{x\,x)G(x,x)G(x,X2) and should enter with the multiplicity 12 
(the number of all pairings of 6 vertices x\,X2,x,x,x,x in which x\ is not connected 
to xi). Let us now check this directly. Indeed, the coefficient of h in (xi, X2) is 

Jjh^)^W el2UGJ)[ ^ = U I dx G te>*M*>*W*>*»)> 

where we applied Wick's theorem to obtain the desired equality. 

In a similar way, the linear term in the expansion of Zjj should correspond to 
the graph with no legs and one vertex of valence four (see Figure [3J| , representing 
J dx G 2 (x,x) (and entering with the multiplicity 3). 

Exercise 4.3 (</> 3 -theory). Let U(<j>) = J dx <fi 3 (x). Find the Feynman rules for 
this theory. Which graphs will contribute to the coefficient of H 2 in the power 
series expansion of the 2-point function {xi,x%)u^ Write down these coefficients 
explicitly. 

4.5. Convergence. Usually the integrals which we get by a perturbative Feynman 
expansion are divergent and ill-defined in many ways. Often one has to renormalize 
(i.e. to find some way to remove divergencies in a unified and consistent manner) 
the theory to improve its behavior. Until recently renormalization was consid- 
ered by mathematicians more like a physical art than a technique; lately Comics 
and Kreimer have done some serious work to explain renormalization in purely 
mathematical terms (see a paper by Kreimer in this volume). 

But even in the best cases, the Green function G{x, y) usually blows up near 
the diagonal x = y, which brings two problems: Firstly, the weights of graphs with 
looped edges, starting and ending at the same point (so-called tadpoles) are ill- 
defined and one has to get rid of them in one or another way. Secondly, all diagonals 
have to be cut out from the spaces over which the integration is performed, so the 
resulting configuration spaces are open and the convergence of all integrals defining 
the weights has to be proved. Mathematically these convergence questions usually 
boil down to the existence of a Fulton-MacPherson-type (see |14| ) compactification 
of configuration spaces, to which the integrand extends. 

There is also a challenging problem to interpret the Feynman diagrams series 
in some classical mathematical terms and to understand the way to produce them 
without a detour to physics and back. In many examples this may be done in terms 
of a homology theory of some grand configuration spaces glued from configuration 
spaces of different graphs along common boundary strata. 

We will see all this on an example of the Chern-Simons theory in the next section. 

5. An example of QFT: Chern-Simons theory 

The Chern-Simons theory has an almost topological character and as such presents 
an interesting object for low-dimensional topologists. For a connection a in a trivial 
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5t/(iV)-bundle over a 3-manifold M one may define (0; see also ^J) the Chern- 
Simons invariant CS(a) as described in Section fo.ll below. It is the action functional 
of the classical Chern-Simons theory and was extensively used in mathematics for 
many years to study properties of 3-manifolds (mostly due to the fact that the 
classical solutions, i.e. the critical points of CS(a), are flat connections). But it 
is the corresponding quantum theory which is of interest for us. Its mathematical 
treatment started only about a decade ago, following Witten's suggestion J2H] that 
it leads to some interesting invariants of links and 3-manifolds, in particular, to 
the Jones polynomial. While Witten's idea was based on the validity of the path 
integral formulation of the quantum Chern-Simons theory, his work catalyzed much 
mathematical activity. By now mathematicians more or less managed to formalize 
the relevant perturbative series and exorcize from them all physical spirit, leaving 
a (surprisingly rich) rigorous mathematical extract. In this section I will describe 
this process in a number of iterations, starting from an intuitive and roughest de- 
scription and slowly increasing the level of rigor and details. Finally, I will try to 
reinterpret these Feynman series in some classical topological terms and formulate 
some corollaries. 

5.1. Chern-Simons theory. Further we will use the following data: 

• A closed orientable 3-manifold M with an oriented framed link L in M. 

• A compact connected Lie group G with an Ad-invariant trace Tr : g — > R 
on the Lie algebra g of G. 

• A principal G-bundle V — > M. 

To simplify the situation, we will additionally assume that G is simply connected, 
since for such groups any principal G-bundle over a manifold M of dimension < 3 
(which is our case) is trivializable, see e.g. jllj . 

The appropriate notions of the Chern-Simons theory, considered as a field theory, 

are as follows. The manifold M plays the role of the space-time manifold X. Denote 

by A the space of G-connections on V and let Q = Aut(P) be the gauge group. 

Fields 4> on M are G-connections on V , i.e. T = A. The Lagrangian is a functional 

L: A-> n 3 (M) defined by 

2 
L(a) = Tr(a A da -\ — a A a A a). 
o 

Remark 5.1. This choice can be motivated as follows. Let = da + a A a be 
the curvature of a. Then Tr(# A 9) is the Chern-Weil 4- form 3 on V, associated 

with Tr; this form is gauge invariant and closed. The Chern-Simons Lagrangian 

2 
CS(a) — Tr(a A H — a A a A a) is an antiderivative of Tr(0 A 6) on V: it is a nice 

exercise to check that d(CS(a)) = Tr(0 A 8). 

The corresponding Chern-Simons action is a function CS : A — > K given by 

If 2 

CS(a) = — / dx Tr(a A da + -a A a A a). 
47r Jm 3 

It is known that the critical points of this action correspond to flat connections and 
(assuming that Tr satisfies a certain integrality property 4 , which holds in particular 



Chern-Weil theory states that the de Rham cohomology class of this form is a certain char- 
acteristic class of V 

Namely that the closed form g— Tr(a A a A a) represents an integral class in H 3 (G, R) 
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for the trace in the fundamental representation of G) it is gauge invariant modulo 
2ttZ. 

The partition function is given by the following path integral: 

(26) Z = [ e ikCS(a) Va. 

J A 

Here the constant k £ N is called level of the theory; its integrality is needed for 
the gauge invariance of Z . 

Now, let L — \J" l =1 Lj, j = 1, . . .m be an oriented framed m-component link in 
M such that each Lj is equipped with a representation Rj of G. Given a connection 
a € A., let holz, (a) be the holonomy 

(27) hol Lj (a) = cxp f a 

of a around Lj. Observables in the Chern-Simons theory are so-called Wilson loops. 
The Wilson loop associated with Lj is the functional 

W(£ J -,i? J )=Tr il ,(hol i ,(a)). 

The m-point correlation function (L) = (Li, L2, . . . L. m ) is denned by 

/in 
e lkcs ^ Y[w{Lj,Rj)Va. 

Since the action is gauge invariant, extrema of the action correspond to points on 
the moduli space of flat connections. Near such a point the action has a quadratic 
term (arising from a A da) and a cubic term (arising from a A a A a). We would 
like to consider a perturbative expansion of this theory. 

5.2. What do we expect. Which Feynman graphs do we expect to appear in the 
perturbative Chern-Simons theory? 

Firstly, a gauge-fixing has to be performed, so the ghosts have to be introduced. 
As a result, we should have two types of edges: the usual non-directed edges (cor- 
responding to the inverse of the quadratic part) and the directed ghost edges. 

Secondly, in addition to the quadratic term the action contains a cubic term, so 
the internal vertices should be trivalent. Also, this time the cubic term is given by 
an antisymmetric tensor instead of a symmetric one, so one should fix a cyclic order 
at each trivalent vertex, with its reversal negating the weight of a graph. Two types 
of edges should lead to two types of internal vertices: usual vertices where three 
usual edges meet, and ghost vertices where one usual edge meets one incoming and 
one outgoing ghost edge. 

Thirdly, note that the situation with legs is somewhat different from our earlier 
considerations. Indeed, the legs (i.e. univalent ends of usual edges) of Feynman 
graphs, instead of being fixed at some points, should be allowed to run over the 
link L, with each link component entering in (L) via its holonomy (|27|l . To reduce 
this to our previous setting, we can use Chen's iterated integrals to expand the 
holonomy in a power series where each term is a polynomial in a. In terms of a 
parametrization Lj : [0, 1] — > R 3 , this expansion can be written explicitly using the 
pullback L*a of a to [0, 1] via Lj-. 
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hol Lj (a) = 1+1 (L*a)(t) + f (L*a)(t 2 ) A {L*a)(h)+ 

Jo<t<i Jo<t 1 <t 2 <i 

(L*a)(t k ) A ■ ■ ■ A (L*a){h) + 



0<ti<-<tfc<l 

where the products are understood in the universal enveloping algebra U(g) of g. 
Thus we should sum over all graphs with any number kj of cyclically ordered legs 
on each Lj, and integrate over the positions 

(x(ti),...,x(t kj )) 6 1-', 0< h < ■■■ < t kj < 1 

of these legs. 

These simple considerations turn out to be quite correct. Of course, one should 
still find an explicit formulas for the weights of such graphs. An explicit deduction 
of the Feynman rules for the perturbative Chern-Simons theory is described in 
details in PI El ■ Let me skip these lengthy calculations and formulate only the 
final results. For simplicity I will consider only an expansion around the trivial 
connection in M = R 3 . 

5.3. Feynman rules. It turns out (see 0) that the weight system W cs of the 
perturbative Chern-Simons theory splits as W cs = W G W, where W G contains all 
the relevant Lie-algebraic data of the theory (but does not depend on the location 
of the vertices of a graph), and W contains only the space-time integration. Since 
the whole construction should work for any Lie algebra, one may encode the anti- 
symmetry and Jacobi relations already on the level of graphs, changing the weight 
Wp of a graph T to a "universal weight" [r], which is an equivalence class of F 
in the vector space over Q generated by abstract (since we do not care about the 
location in K 3 of their vertices) graphs, modulo some simple diagrammatic anti- 
symmetry and Jacobi relations, shown on Figure 0J The same relations hold for 
graphs with either usual or ghost edges, so we may think that the relations include 
the projection making all edges of one type. 



Figure 4. Antisymmetry and Jacobi relations 

The drawing conventions merit some explanation. It is assumed that the graphs 
appearing in the same relation are identical outside the shown fragment. In each 
trivalent vertex we fix a cyclic order of edges meeting there; unless specified other- 
wise, it is assumed to be counter-clockwise. The edges are shown by dashed lines, 
and the link component Lj (fixing the cyclic order of the legs) by a solid line. An 
important consequence of the antisymmetry relation is that for any graph T with a 
tadpole (a looped edge) we have [r] = due to the existence of a "handle twisting" 
automorphism, rotating the looped edge. Thus from the beginning we can restrict 
the class of graphs to graphs without tadpoles. 
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It remains to describe the weight Wt of a graph I\ Roughly speaking, for each 
internal vertex of T we are to perform integration over its position in R 3 (for a ghost 
vertex we should also take a certain derivative acting on the term corresponding 
to the outgoing ghost edge); for each leg we are to perform integration over its 
position in Lj (respecting the cyclic order of legs on the same component). As for 
the edges, we are to assign to each usual and ghost edge inverses of the operator 
curl and of the Laplacian, respectively. 

Somewhat surprisingly (see e.g. \Ej\) two types of edges may be neatly joined 
into one "combined" edge, thus reducing the graphs in question to graphs with only 
one type of edges (and just one type of uni- and trivalent vertices). The weight 
G(x(e),y(e)) of such an edge e with the ends in (x(e),y(e)) £ I 3 x I 3 has a nice 
geometrical meaning: it is given by G(x, y) = lu(x — y), where 

x dx A dx 
uj(x) = — — r: — — h cyclic permutations of (1,2,3) 

2-k\\x\\ 6 

is the uniformly distributed area form on the unit 2-sphere \x\ = 1 in the standard 
coordinates in R 3 . In fact, the usual and the ghost edges (with two possible orien- 
tations) give respectively the (1, 1), (2,0), and (0,2) parts of u>(x — y) in terms of 
its dependence on dx and dy. Abusing notation, I will depict the combined edge 
again by a dashed line. 

Remark 5.2. A simple explanation for an existence of such a simple unified propa- 
gator escapes me. The only explanation which I know is way too complicated: it is 
the existence (see |2]) of the "superformulation" of the gauge-fixed theory, i.e. the 
fact that the connection together with the ghosts may be united in a "superconnec- 
tion" of a supertheory, which leads to an existence of a "superpropagator" , uniting 
the usual and the ghost propagator. I believe that there is a simple explanation, 
probably emanating from the scaling properties and the topological invariance of 
the Chern-Simons theory, by which one should be able to predict that the combined 
propagator should be dilatation- and rotation- invariant. 

Remark 5.3. Note that the weight G(x,x) of a tadpole is not well-defined, so it is 
quite fortunate that we got [r] = for any such graph. 

To sum it up, we are interested in the value 

r 

where |T| is half of the total number of vertices (univalent and trivalent) of T, and 
the weight of T is given by the integral 

(30) Wr= f l[G(x(e),y(e)) 

Jc v e 

over the space Cr of all possible positions of vertices of T, such that all vertices 
remain distinct. Here h = (k + h v )~ 1 , where /i v is the dual Coxeter number of G 
(see [2H1 )■ 

I shall describe in more details the type of graphs which appear in this formula 
and their weights Wt (both the configuration spaces Cr, and the integrand). 
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5.4. Jacobi graphs. Let us start with the graphs. Instead of thinking about 
graphs embedded in R 3 , consider abstract graphs (with just one type of edges), 
such that 

• all vertices have valence one (legs) or three; 

• there are no looped edges; 

• all legs are partitioned into m subsets l\, . . . , l m ; 

• legs of each subset lj are cyclically ordered; 

• each trivalent vertex is equipped with a cyclic order of three half-edges 
meeting there; 

for technical reasons it will be convenient to think that, in addition to the above, 

• all edges are ordered and directed. 

We will further address the last three items simply as an orientation of a graph. 

For such a graph T with a total of 2n (univalent and trivalent) vertices define the 
degree of T by \T\ = n, and denote the set of all such graphs by Zn ■ Set 3 = U3 n ■ The 
ordering and directions of edges of graphs in 3 may be dropped by an application 
of an obvious forgetful map. See Figure |S] for graphs of degree one with m = 2 
and m = 1, and graphs of degree two with m = 1. Both antisymmetry and Jacobi 
relations of Figure 0] preserve the degree of a graph, thus we may consider a vector 
space over Q generated by graphs in 3n modulo forgetful, antisymmetry and Jacobi 
relations. We will call it the space of Jacobi graphs of degree n and denote it by -3„; 
denote also 3 = ®n3n> an d let as before [F] be the class of T E 3 in 3- 



o 







Figure 5. Graphs of degree one and two 

Exercise 5.4. Let m = 2. Write the relations between the equivalence classes of 
degree two graphs shown in Figure Et. What is the dimension of ^2? 

This settles the type of graphs appearing in formula (|29|) : the summation is over 
all graphs in 3, while (L) 6 3[[7i]]- It is somewhat simpler to study separately the 
components of different degrees; define 



(31) (L)n=^ TT ^[T] 

rear, 



^ |Aut(r)| 



5.5. Configuration spaces. Let us deal now with the weights (|30|l of graphs (see 
[51 1181 04] for details). The domain of integration in (|30fl is the configuration space 
CV of embeddings of the set of vertices of T to R 3 , such that the legs of each subset 
lj lie on the corresponding component Lj of the link L in the correct cyclic order. 
It is easy to see that for a graph T with k trivalent vertices and kj legs ending 
on Lj, j = l,...,m we have CV = (M 3 ) fc x n - (S rl x o- kj ~ 1 ) \ A, where cr k is 
a fc-dimensional simplex, and A is the union of all diagonals where two or more 
points coincide. Indeed (forgetting for a moment about coincidences of vertices), 
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each trivalent vertex is free to run over R 3 , while kj legs ending on Lj run over 
S 1 x a l , where S* 1 encodes the position of the first leg, and the following legs are 
encoded by their distance from the previous one. 

Exercise 5.5. Show that the dimension of Cr is twice the number of the edges of 

r. 

Now, an orientation of a graph T determines an orientation of Cr; its idea is in 
fact based on Exercise 15 .51 Let me describe this construction in some local coordi- 
nates. Near each trivalent vertex of T there are three local coordinates (describing 
its movement in R 3 ); assign one of them to each of the three ends edges meeting 
in this vertex using their cyclic order. Near each leg of T there is only one local 
coordinate (describing its movement along the link); assign it to the corresponding 
end of the edge. By now the end of any edge has one coordinate assigned to it. It 
remains to order them using the given ordering of all edges of F and their directions. 
Let us order them as (xi, 2/1,^2,2/2) •• • ,x n ,y n ) where (xi,yi) are the coordinates 
assigned to the beginning and the end of i-th edge. This defines an orientation of 
C r . 

Exercise 5.6. The above construction involves a choice in each trivalent vertex since 
we had only a cyclic order of the edges meeting there, while we used a total order 
of these three edges. Show that a cyclic permutation of the three local coordinates 
used there preserves the orientation of Cr- Also, we used the orientation of T; what 
happens with the orientation of Cr if: 

(1) The cyclic order of three half-edges in one vertex is reversed? 

(2) A pair of edges is transposed in the total ordering of all edges? 

(3) The direction of an edge is reversed? 

5.6. Gauss-type maps of configuration spaces. To understand the integrand 
in l|30|) . consider a directed edge e. Its ends (x,y) represent a point in the square 
R 3 x R 3 with the diagonal A = {(x,y)\ x = y} cut out. This cut square C = 
R 3 x R 3 \ A has the homotopy type of S 2 , with the Gauss map 

<p : (x,y) i-> - - 

lly-^ll 

providing the equivalence. The form uj(y — x) assigned to this edge is nothing more 
than a pullback of the area form o> on S 2 to C via the Gauss map: 

Uj(y — x) = (j)*LU 

Each edge e of a graph T defines an evaluation map ev e : Cr — > C, by erasing 
all vertices of T but for the ends of e. The composition </> e = </> o ev e defines the 
Gauss map corresponding to (the ends of) an edge e. The graph T with an ordering 
ei, . . . , e„ of edges defines the product ^r = II"=i 0e 4 : Cr — ► (S 2 ) n of Gauss maps. 
Finally, the weight Wt is given by integrating the pullback of the volume form 
dvol = A" =1 w on (S 2 ) n to Cr by the product Gauss map </>?: 



(32) W v = / # dvol 

JCt 

Exercise 5.7. Suppose that a graph T has a double edge (i.e., a pair of edges both 
endpoints of which coincide). Show that dim(0r(Cr)) < dim(Cr) — 1. Deduce that 



24 MICHAEL POLYAK 

The following important example shows that at least in some simple cases Wt 
has an interesting topological meaning: 

Example 5.8. Let T = e be a graph with one edge with the ends on two link 
components L, see Figure EK- The configuration space C e = S 1 x S 1 C C is a 
torus. It is mapped to S 2 by the Gauss map <f> — 4> e . The weight W e — J c 4>*uj — 
deg((/)) is in this case just the degree of the map <j>. This fact has many important 
consequences. In particular W e takes only integer values and is preserved if we 
change the uniformly distributed area form u to any other volume form dvol on S 2 
normalized by J„ 2 dvol — 1. It is also preserved if we deform the link by isotopy 
(since then the configuration space changes smoothly and the degree can not jump), 
so is a link invariant. This invariant is easy to identify: J c 4>*uj — lk(Li,L2) is 
the famous Gauss integral formula for the linking number lk(Li, L2) of L\ with L^. 
Thus we get 

Proposition 5.9. Let T — e be a graph with one edge with the ends on two link 
components L. Then the weight W e is the linking number lk(Li, L-i). 






a b 

Figure 6. Signs of crossings and the south pole on S 2 

Exercise 5.10. There is a simple combinatorial way to compute lk(Li,L2) from 
any link diagram: count all crossings where L\ passes over L2, with signs shown in 
FigureEK- Interpret this formula as a calculation of deg(^) by counting (with signs) 
the number of preimages of a certain regular value of 4> (hint: look at Figure Eh) • 
What formula would we get if we counted the preimages of the north pole? 

For other graphs the situation is more complicated. For example, let T = e be 
the graph with one edge, both ends of which end on the same link component, see 
FigureEp- Then the configuration space C e is an open annulus (R 3 ) x S 1 x a 1 \ A = 
S 1 x (0, 1) (torus cut along the diagonal). The Gauss map <j) is badly behaved near 
the diagonal, so the integrand blows up near the diagonal and we can not extend 
it to the closed torus. The integral nevertheless converges; one way to see it is to 
compactify C e , cutting out of it some small neighborhood of the diagonal. This 
makes C e into a closed annulus Cf = S 1 x [e, 1 — e] (thus making the integral 
convergent) and we can recover the initial integral by taking e — * 0. But the 
Gauss integral W e is no more a knot invariant: it may take any real value under 
a knot isotopy. A detailed discussion on this subject may be found in |S]. Why 
does this happen? The reason is that the compactified space C e is not a torus, 
but an annulus, so has a boundary and the degree of the Gauss map is not well- 
defined. When both ends of the edge start to collide together, the direction of the 
vector connecting them (which appears in the Gauss map) tends to the (positive 
or negative) tangent direction to the knot. The image of the unit tangent to the 
knot under the Gauss map is a certain curve 7 on S 2 . One of the boundary circles 
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S 1 x e and S 1 x (1 — e) of Cf is mapped into 7, while the other is mapped into —7, 
and the weight W e is part of the area of S 2 covered by the annulus </>(Cf ) between 
these curves. Unfortunately, 7 may move on S 2 under an isotopy of L, so this area 
may change. 

In this particular case there is a neat way to solve this problem: let L be framed 
(i.e. fix a section of its normal bundle). We may think about the framing as about 
a unit normal vector n{x) in each point a; of a knot. This allows us to slightly 
deform the Gauss map: (f)fX,y) — > (f>(x,y) +en(y). Now both boundary circles of 
the annulus C e map into the same curve on S 2 (why?) and we may glue the annulus 
into the torus so that the map (f> £ extends to it. It makes W e into an invariant of 
framed knots, called the self- linking number (the same result may be obtained by 
slightly pushing L off itself along the framing and considering the linking number 
of the knot with its pushed-off copy) . 

It turns out that for other graphs there are also no divergence problems, so 
all integrals Wt converge, and that a collision of all vertices of a graph to one 
point (so-called anomaly, see ^lEj) is the only source of non-invariance, exactly 
as for W e above. Thus there is a suitable normalization of the expression (|31|) 
for (L) — ^2 n {L) n which gives a link invariant. To avoid a complicated explicit 
description of this normalization, let me formulate this result as follows: 

Theorem 5.11 ([I], [THj, [21 )• Let L = U™ 1 L l be a link. Then (L) depends only 
on the isotopy class of L and on the Gauss integrals W e (Li) of each component 
Li. In particular, an evaluation of (L) at representatives of L for which W e (Li) = 
■ ■ ■ = W e (L m ) — is a link invariant. 

Remark 5.12. It is known that this is a universal invariant of finite type. In partic- 
ular this means that it is stronger than both the Alexander and the Jones polyno- 
mials (it contains the two-variable HOMFLY polynomial) and all other quantum 
invariants. Conjecturally the anomaly vanishes and this invariant coincides with 
the Kontsevich integral, see |18| . 

Example 5.13. Let L be a knot, and take n — 2. There are four graphs of degree 
two, shown in Figure [SJ;. We will denote the first of them X, and the second by 
Y. By Exercise 15. 71 the weight of the third graph vanishes. Also, choose a framing 
of L so that the self-linking is 0; then the contribution of the last graph vanishes 
(another way to achieve the same result is to add to (L) 2 a certain multiple of the 
self-linking number squared); we can set then [X] = [Y]. Thus we will consider 
simply 

v 2 = - 4>* x (lu Auj) + - <j>y(u Auu Auj). 



4 Jcx 3 JCy 

The first integral is 4-dimensional, while the second is 6-dimensional; none of them 
separately is a knot invariant (see |2l)| for a discussion); however, their sum v 2 is 
(see 0, [20])! This invariant is, up to a constant, the second coefficient of the 
Alexander-Conway polynomial. See |20] for its detailed treatment as the degree of 
a Gauss-type map. 

5.7. Degrees of maps. How can we explain the result of Theorem 15. Ill 1 ' We may 
try to repeat the reasoning of Example 15.81 in the general case. Recall that by 
Exercise 15.51 the dimensions of Cr and (S 2 ) n match, so if Cr would be a closed 
manifold, then equation i[52|) would be a formula for a calculation of the degree of 
0r- In other words, if Cr would have a fundamental class, l|52|) would be its pairing 
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with the pullback cfy^dvol. That would be great: we would know that Wt takes only 
integer values, and would be able to compute it in many ways, including a simple 
counting of preimages of any regular value of <fi. 

Unfortunately, the reality is much worse: Cr is an open space, so the degree 
is not well-defined and even the convergence is unclear. To guarantee the conver- 
gence, we should construct a compactification Cr of Cr to which <f>r extends. This 
however will cause new problems: the space Cr will have many boundary strata. 
There are various way to deal with them: we can relativize some of them (i.e. 
consider a relative version of the theory), cap-off some others (i.e. glue to them 
some new auxiliary configuration spaces), or zip them up (gluing a stratum to itself 
by an involution). But in general, some boundary strata will remain; indeed, in 
Example 15. 131 we have seen that none of Wx or Wy separately can be made into 
a knot invariant. The remedy would be to glue together the configuration spaces 
for different graphs in 3n along the common boundary strata. This tedious work 
can be done indeed |18l I24j and (up to a certain anomaly correction) one can in- 
terpret (L) as the degree of a certain map $ n from a grand configuration space 
C n to (S 2 ) n . Too many technicalities are involved to describe this construction in 
necessary details, so I refer the interested reader to ^J |2] and will present only a 
brief sketch of this construction. 

The first problem is that initially the dimensions of Cr for various T £ % n do not 
match. E.g., in Example 15. 131 for n = 2 the spaces Cx and Cy have dimensions 
4 and 6 respectively. This can be fixed by considering a product Cr X (S 2 ) k of 
Cr with enough spheres to make the maps 0r x (id) fc to have the same target 
space (S 2 ) N for all T. Now one should do the gluings. When two endpoints of an 
edge e of a graph T collide, the corresponding boundary stratum of Cr looks like 
Cq x S 2 for G = r/e. Thus we can glue together such strata for all pairs (r, e) 
with isomorphic G = r/e. Some more, so-called hidden, strata remain after these 
main gluings. Fortunately, each of them can be zipped-up (i.e. glued to itself by a 
certain involution) . The only codimcnsion one boundary strata which remains after 
all these gluings are the anomaly strata, where all vertices of a graph r e 3 collide 
together. These problematic anomaly strata can be glued |18] to a new auxiliary 
space. One ends up with a grand configuration space C n endowed with a map 
$ : C n — > (S 2 ) N (glued from the corresponding maps <fir ■ Cr — > (S 2 ) N ). One may 
show that the cohomology H 2N (C n ) of this space projects surjectively to Zn (see ^7] 
for a similar case of 3-manifold invariants). Then (L) n = (L) n + anomaly correction 
can be interpreted as the degree of <&„, or more exactly, the image in 3n of the 
fundamental class [(S' 2 ) Ar ] under the induced composite map 7ro$* : H 2N (S 2 ) N — > 

H (C n ) — > On- 

5.8. Final remarks. There remain many questions: which compactification should 
we take, why do the antisymmetry and Jacobi relations appear in the cohomology 
of the grand configuration space, etc. Each of them is quite lengthy and is out of 
the scope of this note. We refer the interested reader to [Hl^lE]- A mathematical 
treatment of invariants of 3-manifolds arising from the Chern-Simons theory was 
done in [21 0] ; I especially recommend ^7J . While I do not know whether similar 
Feynman series arising in other topological problems always have a reformulation 
in terms of degrees of maps of some grand configuration space, it seems quite plau- 
sible. There are at least some other notable examples, see e.g. ^Hj for a similar 
interpretation of Kontsevich's quantization of Poisson structures. 
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